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O ' Abstract 

' In this article, we study the pseudoscalar bound state 77^/0 (980) (irrespective of the 

, hadro-charmonium and the molecular state) with the QCD sum rules. Considering 

fH ' the SU (3) symmetry of the light flavor quarks and the heavy quark symmetry, we also 

! study the bound states 77^(400 - 1200), ry^" /o(980) and ?7;,"cr(400 - 1200), and make 

' ' ' reasonable predictions for their masses. 

^ ! 
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1 Introduction 

^ : 

In 2007, the Belle collaboration observed two resonant structures in the vr+vr ip' invariant 

ff^ . mass distribution in the cross section for the process e^e~ — t- 'it'^tt~'iP' between threshold 

'i^ [ and \/s = 5.5 GeV using 673 fb~^ of data on and off the T(45) (T'") resonance, one 

2 ; at 4361 ± 9 ± 9 MeV with a width of 74 ± 15 ± 10 MeV, and another at 4664 ± 11 ± 

O ■ 5 MeV with a width of 48 ± 15 ± 3 MeV (they are denoted as the y(4360) and y(4660) 

I respectively) pQ. The structure y(4660) is neither observed in the initial state radiation 

00 ■ (ISR) process e"'"e~ — "lisR'^^'^~ J /i^ [2]) nor in the exclusive cross processes e"'"e~ — >• 

^ '. DD, DD*,D*D*,DD7r, J/'tPD^*W(*^ |3 H El El El • 

, There have been several canonical charmonium interpretations for the y(4660), such 

I as the 5^5"! state [8], the 6^5i state [9], the 5^5i —A^Di mixing state [10], and some exotic 

' interpretations, such as the radial excited state of the ^(| AcAc) + |S^E[!)) jTT], the vector 



^ ■ cscs tetraquark state [T2j, etc. 



A critical information for understanding the structure of the charmonium-like states 
is wether or not the vrvr comes from a resonance. There is some indication that only the 
y(4660) has a well defined intermediate state which is consistent with the scalar meson 
/o(980) in the tttt invariant mass spectra [T3] . 

In Refs.[lH[T5], Voloshin et al argue that the charmonium-like states y(4660), Z(4430), 
y(4260), • • • may be hadro-charmonia. The relatively compact charmonium states {J/ip, 
V'' and Xcj) can be bound inside light hadronic matter, in particular inside higher reso- 
nances made from light quarks and (or) gluons. The charmonium state in such binding 
retains its properties essentially, the bound system (hadro-charmonium, a special molecu- 
lar state) decays into light mesons and the particular charmonium. 

In Ref . p6] . Guo et al assume that the y(4660) is a V''/o(980) bound state (molecular 
state), as the nominal threshold of the tp' — /o(980) system is about 4666 ± 10 MeV [17j, 
the y(4660) decays dominantly via the decay of the scalar meson /o(980), y(4660) — 
ip'f o{980) — )• ip'iTTT, ip'KK, the difficulties in the canonical charmonium interpretation is 
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overcome. Considering the heavy quark spin symmetry, Guo et al predict an r/^/o(980) 
bound state (y(4616)) as the spin-doublet partner of the Y(4660) with a mass of 46161:^ MeV 
and a width of 60 it 30 MeV through the prominent decay mode ij'^tttt [18j . 

In previous work [1^, we studied the mass of the y(4660) as a ijj' fo (980) bound state 
(irrespective of the hadro-charmonium and the molecular state) using the QCD sum rules 
[201 m]. In this article, we extend our previous work to study the r/^/o(980) bound state 
y(4616), furthermore, we take into account the SU{3) symmetry of the light flavor quarks 
and the heavy quark symmetry, study the related hidden charm and hidden bottom states. 
In the QCD sum rules, the operator product expansion is used to expand the time-ordered 
currents into a series of quark and gluon condensates which parameterize the long distance 
properties of the QCD vacuum. Based on the quark-hadron duality, we can obtain copious 
information about the hadronic parameters at the phenomenological side \20\ [2T] . 

The article is arranged as follows: we derive the QCD sum rules for the pseudoscalar 
charmonium-like state y(4616) and the related bound states in section 2; in section 3, 
numerical results and discussions; section 4 is reserved for conclusion. 

2 QCD sum rules for the y(4616) and related bound states 

In the following, we write down the two-point correlation functions Il{p) in the QCD sum 
rules, 



where the Q denotes the heavy quarks c and b. We use the currents J(x) and r]{x) {Q = c) 
to interpolate the bound states r/^/o(980) (predicted in Ref.[l8j) and ry^o"(400 — 1200), 
respectively. 

The hidden charm current c{x)ij5c{x) can interpolate the charmonia r]c, rj'^, , • • • ; 
and the hidden bottom current b{x)i'y5b{x) can interpolate the bottomonia iji,, r(^^ r]'^, ■ ■ ■ . 
We assume that the scalar mesons /o(980) and (t(400— 1200) are the conventional qq states, 
to be more precise, they have large qq components, while in Refs.[T6l I18| the scalar meson 
/o(980) is taken as the KK molecular state. There are hot controversies about their nature, 
for example, the conventional qq states (strongly affected by the nearby thresholds) [22] . 
the tetraquark states, the molecular states [23l[2l]. The currents J{x) and r](x) [Q = c) 
have non- vanishing couplings with the bound states ??c/o(980), r/^/o(980), r/"/o(980), • • • 
and ?7cCt(400 - 1200), r/;.o-(400 - 1200), r?"o"(400 - 1200), • • • , respectively. Considering 
the heavy quark symmetry, there maybe exist some hidden bottom bound states, for 
example, %/o(980), 7y^/o(980), </o(980), <7o(980), %a(400 - 1200), r/^a(400 - 1200), 
r/^V(400 - 1200), r/^'V(400 - 1200), • • • , we study those possibilities with the currents J{x) 
and r][x) {Q = b), and make predictions for their masses. 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J(x) and 'r](x) into the correlation functions Il{p) to 
obtain the hadronic representation [20^ [21] . After isolating the ground state contributions 




(1) 



J(x) 




Q{x)i'y5Q{x)s{x)s{x) , 
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from the pole terms of the Y, we get the following result, 

My — p 

where the pole residues (or couplings) Ay are defined by 

Ay = (0|J/77(0)|y(p)). 



(3) 



(4) 



The contributions from the two-particle and many-particle reducible states are small 



enough to be neglected |25], for example, 



iX: 



Vcfo 



ni: 



Vci 



(p - qf 



m 



fo 



+ ••• , 



(5) 



where the pole residue (or coupling) A^^jp is defined by 

{o\JmvcMp)) = Kfo- 



(6) 



The coupling Xricfo written in terms of the rjc meson decay constant /^^ and the 

coupling Ajq of the scalar meson /o(980) with a tetraquark current. The coupling Ajg 
should be very small as the /o(980) is a light flavor meson, and the two-particle reducible 
contributions can be neglected \25\ I26|. Furthermore, the pseudoscalar charmonia rjc, r]'^, 
T]'^, ■ ■ ■ and the pseudoscalar bottomonia r?^, • • • also have Fock states with addi- 
tional qq components beside the QQ components. The currents J{x) and ri{x) may have 
non-vanishing couplings with the pseudoscalar charmonia and pseudoscalar bottomonia, 
those couplings are supposed to be small, as the main Fock states of the pseudoscalar 
charmonia and pseudoscalar bottomonia are the QQ components, and the pseudoscalar 
charmonia and pseudoscalar bottomonia have much smaller masses than the corresponding 
molecular states Y. 

After performing the standard procedure of the QCD sum rules, we obtain two sum 
rules in the cess and bbss channels respectively: 



A^e" 



dsp{s)e 



(7) 



p[s) = pois) + P(,,)(S) + p^c^^{s){^^) + p^,,y.{s) + P^GGOisMfabcC^G'G') (8) 



vr 



The explicit expressions of the spectral densities po{s), P(ss)is), P{gg){^)^ P{ss)^{^)^ 
P{GGG)i^) presented in the appendix. The sq is the continuum threshold parameter 
and the M'^ is the Borel parameter. We can obtain two sum rules in the ccqq and bbqq 
channels with a simple replacement nis — )• ruq, (ss) — )• (qq) and {sggtrGs) — )• {qgsCrGq). 

We carry out the operator product expansion (OPE) to the vacuum condensates 
adding up to dimension-10. In calculation, we take assumption of vacuum saturation 
for high dimension vacuum condensates, they are always factorized to lower condensates 
with vacuum saturation in the QCD sum rules, factorization works well in the large Nc 
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limit. In this article, we take into account the contributions from the quark conden- 
sates {ss), (ss)^, mixed condensates (sggCrGs), {ss){sgs(TGs), {sgsaGs)'^ and the gluon 
condensates ( "'^'^ ) , {g^fabcG°'G^G'^) (one can see the appendix for the explicit expres- 
sions). The contributions from the quark-gluon condensates { °'"^^ ){ss) , {^^^^^) (sggCrGs) , 
{ggfabcG"'G''G^){ss), ( °°^'^ )(gs)^ are suppressed by large denominators and would not 
play any significant roles. Comparing with the gluon condensate ( "°'^'^ ), the vacuum 
condensates {glfahcG^G^G"), {^^){sgsaGs), {glfabcG''G^G^){ss) are of higher order in 

their contributions are greatly suppressed. In calculation, we observe that the contri- 
butions from the term {g^fabcG°'G^G'^) are less than 0.2%, and can be neglected safely. 

Differentiate the Eq.(7) with respect to then eliminate the pole residue Ay, we 
can obtain the sum rule for the mass of the bound state Y , 

Ja dsp{s)e ^ 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 it 0.01 GeV)^, 
{ss) = {0.8±0.2){qq), {qgsaGq) = ml{qq), (s^.^Gs) = ml{ss), ml = (0.8 ± 0.2) GeV^, 
(^) = (0.33GeV)^ {g^sfabcG'^G'G'^) = 0.045 GeV^, m„ = « 0, = (0.14 ± 
0.01) GeV, rric = (1.35±0.10) GeV and ruh = (4.8±0.1) GeV at the energy scale p = I GeV 
[201 EH [27] . 

In the conventional QCD sum rules [20','21|, there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter and 
threshold parameter sq. We impose the two criteria on the pseudoscalar charmonium-like 
states Y to choose the Borel parameter and threshold parameter sq. 

We take the threshold parameter of the pseudoscalar bound state 7/^/o(980) as s^j = 
(4.62 + 0.5)^ GeV^ w 26 GeV^ tentatively to take into account possible contribution from 
the ground state, where the energy gap between the ground state and the first radial 
excited state is chosen to be 0.5 GeV. Taking into account the SU{3) symmetry of the light 
fiavor quarks, we expect the threshold parameter s^^ (for the bound state 7/^(t(400 — 1200)) 
is slightly smaller than the s%. Furthermore, we take into account the mass difference 
between the c and b quarks, the threshold parameters in the hidden bottom channels are 
tentatively taken as Sgg = 142 GeV^ and s% = 144 GeV^. In this article, we use those 
value as a guide to determine the threshold parameters so with the QCD sum rules. 

The contributions from the high dimension vacuum condensates in the operator prod- 
uct expansion are shown in Figs. 1-2, where (and thereafter) we use the (qq) to denote the 
quark condensates {qq), {ss) and the {qgsCrGq) to denote the mixed condensates {qgsCrGq), 
{sggCrGs). From the figures, we can see that the contributions from the high dimension con- 
densates change quickly with variation of the Borel parameter at the values < 2.7 GeV^ 
and < 7.6 GeV^ in the hidden charm and hidden bottom channels respectively, such 
an unstable behavior cannot lead to stable sum rules, our numerical results confirm this 
conjecture, see Fig.4. 

At the values > 2.7 GeV^ and so > 26GeV^, the contributions from the {qq)'^ + 
{qq){qgs(^Gq) term are less than 19.5% in the cess channel, the corresponding contributions 
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are less than 35.5% in the ccqq channel; the contributions from the vacuum condensate of 
the highest dimension (qgsaGq)'^ are less than 7% in all the hidden charm channels, we 
expect the operator product expansion is convergent in the hidden charm channels. At the 
values > 7.6 GeV^ and so > 144 GeV^, the contributions from the {qq)'^ + {qq){qgs(yGq) 
term are less than 7% in the hbss channel, the corresponding contributions are less than 
18% in the bhqq channel; the contributions from the vacuum condensate of the highest 
dimension {qgsuGq)'^ are less than (or equal) 6% in all the hidden bottom channels, we 
expect the operator product expansion is convergent in the hidden bottom channels. 

The contributions from the gluon condensate (^^^^) are less than (or equal) 37% 
(26%) in the cess (ccqq) channel at the values > 2.7 GeV^ and sq > 26 GeV^; while the 
contributions are less than 19.5% (16.5%) in the bbss (bbqq) channel at the values > 
7.6 GeV^ and sq > 144 GeV^. The contributions from the high dimension condensates 
^ OsGG ^ [(gg) + (qgsCrGq) + (qq)"^] are small enough and neglected safely. 

In the QCD sum rules for the tetraquark states (irrespective of the molecule type and 
the diquark-antidiquark type), the contributions from the gluon condensate are suppressed 
by large denominators and would not play any significant roles for the light tetraquark 
states |28| I29j. the heavy tetraquark state [30] and the heavy molecular state [3T]. In the 
present case, the contributions from the gluon condensate ( "°^'^ ) are rather large, just like 
in the sum rules for the y(4660) |19| . If we take a simple replacement s{x)s{x) — t- (ss) and 
[u{x)u{x) + d{x)d{x)~\ — )• 2{qq) in the interpolating currents J(x) and ??(a;), we can obtain 
the standard pseudoscalar heavy quark current Q{x)i^^Q{x)^ where the gluon condensate 
(^asGG_\^ plays an important rule in the QCD sum rules [20]. 

In this article, we take the uniform Borel parameter M^^^, i.e. M^^„ > 2.7 GeV^ and 
^min — '^■^ GeV^ in the hidden charm and hidden bottom channels, respectively. 

In Fig. 3, we show the contributions from the pole terms with variation of the Borel pa- 
rameters and the threshold parameters sq. If the pole dominance criterion is satisfied, 
the threshold parameter sq increases with the Borel parameter monotonously. From 
Fig.3-A, we can see that the pole dominance criterion cannot be satisfied at the values 
So < 25 GeV^ and > 2.7 GeV^ in the cess channel, the threshold parameter so has to 
be pushed to larger value. 

The pole contributions are larger than 48% at the values < 3.1 GeV^ and so > 
25GeV2, 26GeV2 in the ccqq, ecss channels respectively; and larger than 50% at the values 
< 8.2 GeV^, sq > 142 GeV^ 144 GeV^ in the bbqq and bbss channels respectively. 
Again we take the uniform Borel parameter M^^^,, i.e. M^^^ < 3.1 GeV^ and M^^^. < 
8.2 GeV^ in the hidden charm and hidden bottom channels, respectively. 

In this article, the threshold parameters are taken as sq = (26±l)GeV2, (27±l)GeV2, 
(144 ± 2) GeV^ and (146 ± 2) GeV^ in the ccqq, cess, bbqq and bbss channels, respectively; 
the Borel parameters are taken as = (2.7-3.1) GeV^ and (7.6-8.2) GeV^ in the hidden 
charm and hidden bottom channels, respectively. In those regions, the pole contributions 
are about (48 - 72)%, (49 - 72)%, (50 - 66)% and (51 - 66)% in the ecss, ccqq, bbss 
and bbqq channels, respectively; the two criteria of the QCD sum rules are fully satisfied 
[201 EI]. 

The Borel windows M"^^^ — M^-^ change with variations of the threshold parameters 
So, see Fig. 3. In this article, the Borel windows are taken as 0.4 GeV^ and 0.6 GeV^ 
in the hidden charm and hidden bottom channels respectively, they are small enough. 
Furthermore, we take uniform Borel windows and smear the dependence on the threshold 
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parameters sq- If we take larger threshold parameters, the Borel windows are larger 
and the resulting masses are larger, see Fig. 4. In this article, we intend to calculate the 
possibly lowest masses which are supposed to be the ground state masses by imposing the 
two criteria of the QCD sum rules. 

In Fig. 4, we plot the bound state masses My with variation of the Borel parameters 
and the threshold parameters. The hidden charm current c{x)i'y^c{x) can interpolate the 
charmonia r/c, rj'^, rj'^-, and the hidden bottom current h{x)i'^^h[x) can interpolate the 
bottomonia ry^, 7/[,', ••• [T7]. The currents J{x) have non- vanishing couplings with 
the bound states %/o(980), ?7c/o(980), r]'^Wm), • • • and ??fe/o(980), r?^/o(980), r?^7o(980), 
r/^"/o(980), • • • , respectively. The mass of the r/^ listed in the Particle Data Group is 
= 9388.9l2;3 =t 2.7 MeV, while the r/^, r/^', r/^", • • • are not observed yet [17]. In the 
constituent quark models, the mass splitting between the spin-singlet and spin-triplet are 
proportional to ^fj^:, in the heavy quark limit, the r]b and T degenerate. The constituent 
quark mass Mf, is large enough, Mx = (9460.30 ± 0.26) MeV, the energy gap between the 
r/f, and T is about 71.4 MeV, the energy gaps between the radial excited states are even 
smaller. In this article, we assume the masses of the ry^,, r/^', • • • are slightly smaller 
than ones of the T, T', T", • • • respectively. 

From Figs.3-A,3-C,4-A,4-C, we can see that the QCD sum rules support existence of 
the ry^/o(980) and ?7^"/o(980) bound states, the nominal thresholds of the rjc — /o(980) and 
T" — /o(980) systems are too low, and we cannot reproduce the 77c/o(980) and rj'i^f q(980) 
bound states. Our predictions for the masses the Y(4660) [T9] and y(4616) support the 
conjecture of Voloshin et al, i.e. a formation of hadro-charmonium is favored for higher 
charmonium resonances ip' and Xcj as compared to the lowest states J/ip and r]c [14J. 

In this article, we intend to prove that the r/^/o(980) bound state can be reproduced by 
the QCD sum rules, the pseudoscalar charmonium-like state y(4616) predicted in Ref.|18| 
maybe exist. 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the masses and pole resides of the pseudoscalar bound states Y, which are shown 
in Figs. 5-6 and Tables 1-2. In this article, we calculate the uncertainties 5 with the formula 

E(^)'u«=^>(^^-^^)'' (10) 

where the / denote the hadron mass My and the pole residue Ay, the Xi denote the 
input QCD parameters rric, nib, (qq), {ss), • • • , and the threshold parameter sq and Borel 
parameter M^. As the partial derivatives ^ are difficult to carry out analytically, we take 

the approximation ^ (xi—Xi)"^ w [f{xi ± Axi) — /(xj)]^ in the numerical calculations. 

In table 1, we also present the nominal thresholds of the t/^ — /o(980), 77^ — cj(400 — 1200), 
T'" - /o(980) and T" - cr(400 - 1200) systems. From the table, we can see that there 
maybe exist a bound state r7^/o(980) as the partner of the y(4660). The predicted mass of 
the bound state r/^cr(400 — 1200) is about (4.56 ± 0.21) GeV, while the nominal threshold 
of the r/^ - a(400 - 1200) system is about (4.037 - 4.837) GeV. There maybe exist such 
a bound state. The bound states r/^/o(980) and r/^cr(400 — 1200) can be produced in the 
exclusive decays of the B meson through b — )• ccs, ccq at the quark level. 

In the bbss channel, the numerical result My = 11.42 it 0.21 GeV indicates that there 
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bound states 


My (GeV) 


A/,,,/T- + My„/, (GeV) 


cess 


4.68 ±0.29 


4.617 


ccqq 


4.56 ±0.21 


4.037- 4.837 


bbss 


11.42 ±0.21 


11.559 


bbqq 


11.36 ±0.18 


10.979- 11.779 



Table 1: The masses of the pseudoscalar bound states. 



bound states 


Xy (10-^GeV^) 


cess 


3.63 ±1.80 


ccqq 


3.41 ± 1.37 


bbss 


20.7 ±8.0 


bbqq 


20.4 ±6.5 



Table 2: The pole residues of the pseudoscalar bound states. 

maybe exist a bound state T"'/o(980), which is consistent with the nominal threshold 
Mr'" + Mf^ = 11.559 GeV, while the nominal thresholds Mr + Mj^ = 10.44 GeV, Mr' ± 
M/o = 11.00 GeV, Mr" ± Mj^ = 11.335 GeV are too low. The scalar meson c7(400 - 1200) 
is rather broad with the Breit-Wigner mass formula (400 — 1200) — i(250 — 500) |17| . 
Considering the SU{3) symmetry of the light flavor quarks, we can obtain the conclusion 
tentatively that there maybe exist the bound states 77^ct(400 — 1200) and ?7"'cj(400 — 1200) 
which lie in the regions (4.037 - 4.837) GeV and (10.979 - 11.779) GeV, respectively. As 
the energy gaps between the T's are rather small and the scalar meson cr(400 — 1200) is 
broad enough, there maybe exist the 7ybCT(400- 1200), r/[,cT(400- 1200) and 77^'o-(400- 1200) 
bound states. We cannot draw decisive conclusion with the QCD sum rules alone. 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. In 
proton-proton collisions at ^/s = 14TeV, the bb cross section is expected to be ~ 500 fib 
producing 10^^ bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 X 1032(.jj^-2gg(,-i |32j_ Ti^g pseudoscalar bound states r/[,'7o(980) and ry^"cj(400 - 1200) 
predicted in the present work may be observed at the LHCb, if they exist indeed. We 
can search for those bound states in the r/feTrvr, J^^vrvr, r/^Vvr, rf'^'mT, r]hKK, rj'i^KK, tj'j^KK, 
Tj'i^'KK, ■ ■ ■ invariant mass distributions. 

4 Conclusion 

In this article, we study the pseudoscalar bound state 77^/0 (980) (irrespective of the hadro- 
charmonium and the molecular state) with the QCD sum rules, the numerical result 
My = 4.68 ± 0.29 GeV is consistent with the value 4616;g MeV predicted by Guo et al. 
Considering the SU{3) symmetry of the light flavor quarks and the heavy quark symmetry, 
we also study the bound states rj'^a{AOO - 1200), rj'l'f o{980) and r]'^'a{400 - 1200) with the 
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QCD sum rules, and make reasonable predictions for their masses. Our predictions depend 
heavily on the two criteria (pole dominance and convergence of the operator product 
expansion) of the QCD sum rules. We can search for those bound states at the LHCb, 
the KEK-B or the Fermi-lab Tevatron. 

Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 

P^^'^ = daj' '^d0a0{l-a-0)\s-mlf{3s + ml), (11) 



9ms{ss) rf , Z"^-" ,^ ^ , ^2 N 
P{ss){s) = g^^4 J da d0aPs{s-mQ) 



THsisgsaGs) 
327r4 



daa{l — a){3s — mg) , 



(12) 



P{ss)^{s) = J daa{l - a){3s - uiq) 



167r2 

{ss){sgs(TGs) 
32^r2 



daa{l — a) 



6+(4. + ^l<5(s-mQ) 



^ m^Q{ss){sgsaGs) 



327r2 



J ai 



1 + 



M2 



3m|(s5sCrGs)2 



s 

6M6 



7687r2M6 



/ das 6{s — mq) , 



(13) 



P{GG) {s) 



2567r4 



da 



ft 



dPaPs{s — rrir 



5127r4 

2 



/ da d/3(l 

J Jft 



a - (3fs{s - mo) 



10247r4 

^4 



m7 



+ 



10247r4 
10247r4 



da / d/3 

-'ft 

da / d/3 

-'ft 

dp 



/ c/a / 



a j3 
1 1 
1 1 



(1 - a - /3)'(3s - 2m|) 
{l-a-pf 

(1 _ a - /3)2(s - m^) , (14) 
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24 2.6 2.8 3.0 3.2 3.4 3.6 3.8 4.0 42 44 

M" [GeV"] 




M' [GeV'] 




Figure 1 : The contributions from different terms witli variation of the Borel parameter 
in the operator product expansion. The A, B and C correspond to the contributions from 
the (qgsCrGq)'^ term, the (qq)'^ +{qq){qgs(TGq) term and the ( "°^'^ ) term, respectively. 
The (I) and (II) denote the cess and ccqq channels, respectively. The notations a, /3, 
7, A, /9 and r correspond to the threshold parameters sq = 23GeV^, 24GeV^, 25GeV^, 
26GeV2, 27GeV2 and 28GeV2, respectively. 
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Figure 2: The contributions from different terms with variation of the Borel parameter 
in the operator product expansion. The A, B and C correspond to the contributions 
from the {qggCrGq)'^ term, the {qq)"^ +{Q<l){Q9sO'Gq) term and the {^^^^) term, respectively. 
The (I) and (II) denote the bbss and bbqq channels, respectively. The notations a, /3, 7, 
A, p and r correspond to the threshold parameters sq = 138 GeV^, 140 GeV^, 142 GeV^, 
144 GeV^, 146 GeV^ and 148 GeV^, respectively. 
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[GeVl 



24 2.6 2.8 3.0 32 3.4 3.6 3.8 4.0 4.2 4.4 



[GeV] 




Figure 3: The contributions from the pole terms with variation of the Borel parameter 
M^. The A, B, C, and D denote the cess, ccqq, bbss and bbqq channels, respectively. In 
the hidden charm channels, the notations a, /3, 7, A, p and r correspond to the threshold 
parameters sq = 23GeV^ lAGeV"^, 25 GeV^, 26GeV^ 27GeV2 and 28 GeV^ respectively 
; while in the hidden bottom channels they correspond to the threshold parameters sq = 
138 GeV^, 140 GeV^, 142 GeV^, 144 GeV^ 146 GeV^ and 148 GeV^, respectively. 
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Figure 4: The masses of the pseudoscalar bound states with variation of the Borel pa- 
rameter and threshold parameter sq. The A, B, C, and D denote the cess, ccqq, bbss, 
and bbqq channels, respectively. In the hidden charm channels, the notations a, /3, 7, A, p 
and T correspond to the threshold parameters sq = 23GeV^, 24GeV^, 25GeV^, 26GeV^, 
27GeV^ and 28 GeV'^, respectively ; while in the hidden bottom channels they correspond 
to the threshold parameters sq = 138 GeV^, 140 GeV^ 142 GeV^, 144 GeV^ 146 GeV^ and 
148 GeV^, respectively. The ^ and denote the rjc — /o(980) and rj'^ — /o(980) thresholds 
respectively in the cess channel, while in the bbss channel they correspond to T" — /o(980) 
and T'" — /o(980) thresholds respectively. 
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